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In this work, we have consider Kaluza-Klein Cosmology for anisotropic universe where the universe 
is filled with variable modified chaplygin gas (VMCG). Here we find normal scalar field (f> and the 
self interacting potential V((f>) to describe the VMCG Cosmology. Also we graphically analyzed 
the geometrical parameters named statefinder parameters in anisotropic Kaluza-Klein model. Next, 
we consider a Kaluza-Klein model of interacting VMCG with dark matter in the Einstein gravity 
Ch , framework. Here we construct the three dimensional autonomous dynamical system of equations for 

this interacting model with the assumption that the dark energy and the dark matter are interact 
between them and for that we also choose the interaction term. We convert that interaction terms 
to its dimensionless form and perform stability analysis and solve them numerically. We obtain a 
stable scaling solution of the equations in Kaluza-Klein model and graphically represent solutions. 
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PACS numbers: 04.20.Fy, 98.80.Cq, 98.80.Es 



D ' I. INTRODUCTION 

wo: 



The recent observations of the luminosity of type la supernovae [1,2] indicate that the universe is currently 
undergoing an accelerated expansion. The theory of Dark energy having negative pressure is the main respon- 
sible for this scenario. From recent cosmological observations including supernova data [3] and measurements 
of cosmic microwave background radiation (CMBR) [4] it is evident that our present universe is made up 
I of about 2/3 dark energy and about 1/3 dark matter of the whole universe. There are several interesting 
mechanisms such as higher dimensional phenomena [5], Loop Quantum Cosmology (LQC) [6], modified gravity 
[7], Brans-Dickc theory [8], branc- world model [9] and many others which are extremely applicable in recent 
. research to explain the present scenario of the universe. 

, A fundamental theory on higher dimensional model was introduced by Kaluza-Klein [10, 11], where they 
OO ■ considered an extra dimension with FRW metric to unify Maxwell's theory of electromagnetism and Einstein's 
| gravity. After extensive research, scientists are realized that higher dimensional Cosmology may be more useful 
to understand the interaction of particles. Also to study the universe, as our space-time is explicitly four 
dimensional in nature, the 'hidden' dimensions must be related to the dark matter and dark energy. In modern 
physics Kaluza-Klein theory showed its great importance in string theory [12], in supcrgravity [13]and in 
superstring theories [14]. Li et al [15] have considered the inflation in Kaluza-Klein theory. In 1997, Overduin 
and Wesson [16] represented a review of Kaluza-Klein theory with higher dimensional unified theories. String 
cloud and domain walls with quark matter in N-dimensional Kaluza-Klein Cosmological model was presented 
m | by Adhav ct al [17]. Recently some authors [18-24] studied Kaluza-Klein cosmological models with different 

■ dark energies and dark matters. 

_ _ i 

Recently dynamical system theory become a very important tool to study cosmology and astrophysics in 
the framework of general relativity. Some authors [26-30] work on qualitative analysis in dynamical system of 
FRW Cosmology in Brans-Dicke gravity, where they showed that the field equations of FRW may be reduced 
to two dimensional dynamical system whether in presence of scalar potential it could be reduced in three 
dimensional. 

In this work, we have considered Kaluza-Klein Cosmology for Anisotropic Universe where the Universe 
is filled with variable modified Chaplygin gas (VMCG). Here we find normal scalar field 4> and the self 
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interacting potential V((f>) to describe the VMCG Cosmology. Also we graphically analyzed the geometrical 
parameters named statefinder parameters in anisotropic Kaluza-Klein model. In the later section, we consider 
a Kaluza-Klein model of interacting VMCG with dark matter in the Einstein gravity framework. Here we 
construct the three dimensional autonomous dynamical system of equations for this interacting model with the 
assumption that the dark energy and the dark matter are interact between them and for that we also choose 
the interaction term. We convert that interaction terms to its dimensionless form and perform stability analysis 
and solve them numerically. We obtain a stable scaling solution of the equations in Kaluza-Klein model and 
graphically represent solutions. 



II. BASIC EQUATIONS AND SOLUTIONS 

The Kaluza-Klein type metric is given by [25] 

ds 2 = dt 2 - a 2 (dx 2 + dy 2 + dz 2 ) - b 2 dip 2 (1) 

where a and b are function of cosmic time t only. Here we are dealing only with an anisotropic fluid whose 
energy momentum tensor is in the following form: 

T^ = diag[T§,Tl,Ti,Ti] 

We parameterize it as follows: 

T" = diag[p, -p x , -p y , -p z , -p^] 
= diag[l, -u) x , -0J y , -ui z , -0J^]p 

= diag[l, — cj, — cj, — cj, —cu]p (2) 

where p is the energy density of the fluid; p x ,Py,Pz and p^ are the pressures and u> x ,u) y ,u) z and uj^ are the 
directional equation of state (EoS) parameters of the fluid. 

The Einstein field equations, in natural limits (&ttG = c = 1) are 

where g^u^u 1 ' — = (1,0,0,0,0) is the velocity vector; is the Ricci tensor; R is the Ricci scalar, 

T M „ is the energy momentum tensor. In a co-moving coordinate system, Einstein's filled equation (3), for the 
anisotropic Kaluza-Klein spacctime (1), with eq.(2) (for fiat universe i.e., for k — 0) yield 

„ ab „ a 2 

3- + 3-=p 4 
ab a z 

'a b a 2 n ab 
2- + - + - + 2- = -p (5) 
aba- ab 

3^ + 3^ = -p (6) 
a a 1 

where the overhead dot(') denote derivative with respect to the cosmic time t. The energy conservation equation 
in Kaluza-Klein Cosmological Model is as 

p + 4H(p + p) = (7) 
In the case of pure Chaplygin gas obeys an equation of state like [31, 32] 



3 



where p and p are respectively the pressure and energy density and B is a positive constant. Subsequently in 
the case of generalized Chaplygin gas [33-35] obeys an equation of state 



B 



with < a < 1 



and in the case of modified Chaplygin gas [36, 37] obeys an equation of state 



V = Ap- 



B 



with < a < 1 



(9) 



(10) 



where A, B arc positive constants. This equation of state shows that radiation era (when A = 1/3) at one 
extreme (when the scale factor a(t) is vanishingly small) while a AC DM model at other extreme (when the 
scale factor a(t) is infinitely large). 

Guo and Zhang 2005a [38] first proposed variable Chaplygin gas model with equation of state (8), where B 
is a positive function of the cosmological scale factor 'a' i.e., B = B(a). This assumption is reasonable since 
B(a) is related to the scalar potential if we take the Chaplygin gas as a Born- Infeld scalar field Bento et al. 
2003 [39] . There are some works on variable Chaplygin gas model like as Sethi et al. 2005 [40] , Guo and Zhang 
20056 [41]. In 2007, Debnath [42] introduced VMCG with equation of state (10) where B is a positive function 
of the cosmological scale factor 'a' (i.e., B = B(a)) as 



p = Ap 



B(a) 



with < a < 1 



(11) 



where A is a positive constant. Now, we have considered that the anisotropic universe filled with VMCG obeys 
equation of state (10) where B is a positive function of the cosmological scale factors a and b (i.e., B = B(a, b)) 



A B(a,b) o < q, < 1 

Now, for simplicity, assume B(a, b) is of the form 

B(a,b) =B (a 3 b)~ n 
where Bq > and n are constants. From (7) using (12) and (13) we have 



(12) 



(13) 



P 



B (a+1) 



C 



((A+l)(a + l)-n) (a 3 b) n ( a 3 6 )(a+i)(A+i) 



(14) 



where C > is an arbitrary integration constant and (l + A)(l+a) > n, for positivity of first term. Here n must 
be positive, because otherwise, a — > oo implies p — >• oo, which is not the case for expanding Universe. Now, 
we consider the energy density and pressure corresponding to normal scalar field <j) having a self-interacting 
potential V(<p) as follows: 



and 



p* = \i 2 + v(-P) = p = 



B (a + l) 



C 



({A+ l)(a+ 1) -n) (a 3 b) n ( a 3&)(a+i)(A+i) 
B {a + 1) 1 



A 



C 



((A+l){a+l)-n) (a 3 6) r 
B (a 3 b)- n 



( a 3fr)(a+l)(A+l) 



gofa+1) 



C 



((A+X)(a+l)-n) (a 3 i>)™ 1 (a 3 6)C<*+i)(A+i) 

Hence for flat universe (i.e., k = 0), we have 

B {a + l)(a 3 b)~ n 



(A+l) 



l + A-n + a + Aa 



+ C(a 3 by 



(X+l)(a+l) 



(15) 



(16) 
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Fig.3 



Fig. 1 and Fig. 2 show the variations of V against a and b for A — | and A = 1 respectively with 
B = 1, C = 2, n = 3, a = 0.8 and Fig. 3 shows the variation of <j> against t for 
A = |,B = 2, C = 1,ti = 1.1, d = 1, q = 0.5 in the case of the anisotropic universe filled with Variable Modified 

Chaplygin Gas with Kaluza-Klein metric. 



-B (a 3 b) 



3 r\ —n 



3u\-n 



B Q (a + l)(a 3 b) 
1 + A-n + a + Aa 



+ C(a 



dt 



(17) 



and 



V 



-(A- I) 



B (a + l)(a 3 b)- n 
1 + A- n + a + Aa 



C a 



-B (a 3 



B (a + l){a 3 b)- n 
1 + A- n + a + Aa 



+ C (a 3 b) 



-(A+l)(a+l) 



(18) 



From equation (17), we see that the integration is very complicated. So the numerical investigations are 
needed to get the physical interpretation of the scalar field and the corresponding potential. Fig. 1 and Fig. 2 
show the variations of V against a and b for A = | and A = 1 respectively with B = 1, C = 2, n = 3, a = 0.8 
and Fig. 3 shows the variation of <j> against t for A = |, B = 2, C = 1, n = 1.1, d = 1, a = 0.5 in the case of the 
anisotropic universe filled with Variable Modified Chaplygin Gas with Kaluza-Klein metric. From the figures, we 
see that <fi increases due to evolution of the universe. Also the potential V decreases with the scale factors a and b. 

In the paper Gorini et al. 2003 [33]; Alam et al. 2003 [34], the fiat Friedmann model filled with Chaplygin 
fluid has been analyzed in terms of the recently proposed statefinder parameters Sahni et al. 2003 [43]. The 
statefinder diagnostic along with future SNAP observations may perhaps be used to discriminate between 
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different dark energy models. The above statefinder diagnostic pair for higher dimensional anisotropic cosmology 
are constructed from the scale factors a and b as follows: 



H H , r- 1 

r=1 + ''jP + JP md S= 3(^i) 



(19) 



where q is the deceleration parameter defined by q = — 1 — jp- and is the Hubble parameter. Since this 
parameters are dimensionless so they allow us to characterize the properties of dark energy in a model indepen- 
dently. Finally we graphically analyzed geometrical parameters r and s with respect to a and b in the case of 
the anisotropic universe filled with Variable Modified Chaplygin Gas. 
For Kaluza-Klein model with flat universe (i.e., k = 0) we have 



4if + 16H 2 = -(p-p) 



(20) 



Therefore 



Hence 



r = 21 



and 



3H 2 



3 - 



P-P 
3H 2 



7 , +S A a P 

- + (a + l)A- — 
4 p 



4 " / . \ 



(21) 



(22) 



3(f 



15 + H^ [ A + {a + l)A ~p ) + 1P {P ~ Af>) 



(23) 



where p and p are given by (15) and (16). Using field equations, r and s can be expressed in terms of scale 
factors a and b in very complicated form, so we have not presented the explicit expressions here. Fig. 4 and 
Fig. 5 show the variations of r against a and b and the variations of s against a and b and Fig. 6 shows the 
variation of r against s for A = 1,B = 1, C = 2, n = 3, a = 0.8 and H = 1 respectively in the case of the 
anisotropic universe filled with Variable Modified Chaplygin Gas. From the figures, we conclude that r increases 
and s decreases as a and b increase. From figure 6, we see that s increases from some positive value to +00 as 
r decreases from finite positive range. Also after certain stage, s increases from —00 to as r decreases from 
some positive to negative value. 



III. DYNAMICAL SYSTEM ANALYSIS IN KALUZA-KLEIN MODEL WITH 
INTERACTING VMCG AND DARK MATTER 

We consider the universe filled with VMCG as dark energy and dark matter interacting through an interaction 
term. Due to interaction between the two components, the energy conservation would not hold for the individual 
components, therefore the above conservation equation (7) will break into two non-conserving equations: 

Pvmcg ~l~ 4.H (pvmcg ~t~ Pvmcg) Q (^^) 

and 

p dm + 4H(p dm + p dm ) = Q (25) 



Here Q is the interaction term which have dimension of density multiplied by Hubble parameter H . For this 
scenario the suitable choice of Q be Q = AbH p. b is the coupling parameter denoting the transfer strength, 
P = Pvmcg + Pdm and p = p vmC g + Pdm {pdm is very small quantity, somewhere dark matter assumed as 
pressureless quantity) are the total cosmic energy density and pressure respectively which satisfy the above 
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Fig.5 



Fig.6 

Fig. 4 and Fig. 5 show the variations of r against a and b and the variations of s against a and b for 
A — 1, B = 1, C — 2, n = 3, a — 0.8 and H = 1 respectively and Fig. 6 shows the variation of r against s for 
A = 1, B = 1, C = 2, n = 3, a = 0.8 and H — 1 respectively in the case of the anisotropic universe filled with Variable 

Modified Chaplygin Gas. 



conservation equation (7). The subscripts vmcg and dm denote the VMCG and dark matter respectively. 
The dark matter equation of state is assumed to be linear EoS 

Pdm = UdmPdm (26) 

To analyze the dynamical system, we convert the physical parameter into dimensionless form as follows: 



i / 3 7 \ Pvmcn Pvjncq Pd 

\n(a a b), u = -ttt^-, v = . , w = 



AH 2 ' 



AH 2 



AH' 2 



(27) 



The equation of state of the VMCG can be expressed as 

Pv7UC( 



Pvr, 



(28) 



Now we can cast the evolution equations in the following autonomous system of u, v and w in the form: 



du 2 

— = U — V — CU — CW — -U{U — V + (1 — UJdmjW) 

dx 3 



(29) 
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dv v 2 

— = ((1 + a)A — a—)(—u — v — c(u + w)) + n(Au — v) — —v(u — v + (1 — u> dm )w) + 2v (30) 
dx u 3 



dw 
dx 



(1 + UJ dm )w + C{u +W)— -W(ll - V + (1 - LOdm)w) + 2w 



(31) 



For the mathematical simplicity, we work out a = 1 and A — i only. The critical points of the above system 

du dv_ dw 

dx dx dx 



are obtained by putting 4 s = 4^ = ^ = which yield 



3 3c 



2 n- 2(1 +uJdm) 

3(n 2 - 2n(2 + uj dm ) + 4(1 + uj dm + cuj dm )) 
4(n - 2(1 + uj dm )) 

3c 



^crit 



2 — n + 2a;^r] 



(32) 
(33) 
(34) 



For the stability of the dynamical system about the critical point, we linearize the governing equation around 
the critical point and arrived at 



(3 - 3c - 4m + 2v + 2w{-\ + uj dm )) 



6u + 


2u 


Sv + 


1 


- crit 


3 


crit 



{-3c + 2u(-l + cu dm )) 



Sw 

it 

(35) 



{ dv 

dx 



3Au 2 (l + c-n + a + ca)+ v(2u 2 + 3(v + cw)a) 



3u 2 



Su 



2u Av 
2-n- — + — + {! + c)a 



(2v + cw)a . .„ . 2 

^ '- A(l + a) + -w(-l + Lu dm ) 

U 3 



Sv 



. /- \ 1 , „ 3ca 

-Ac(l + a) + -w(-2 + + 2cj dm 

3 it 



5w (36) 





2w~ 




~2w~ 


Sv + 






5u + 




\~dx) ~ 




crit 


3 


crit 



-(3c - 2u + 2v + (-3 + 4tt»)(-l + w dm )) 



The corresponding eigenvalues of the coefficient matrix of the above equations are 

n — 4 



5w 



and 



^2,3 — 



Ai 



44 - 40n + 9n 2 + 4c(3n - 8) + 56cj dm - 24nw dm + 12w 



12(2c + n-2(l+w dm )) 



(37) 



(38) 



± 



y/(8c(3u dm - 1) + (6u dm + 3n - 10)(n - 2(u dm + l)))(-48c 2 + 8c(9uw - 3n + 5) + {6tj dm + 3n - 10)(?t - 2(u dm + 1))) 

12(2c+n-2(l+w dm )) 



(39) 

If the real parts of the above eigenvalues are negative at critical points then critical points are stable and 
are stationary attractor solutions; otherwise unstable and thus oscillatory. Therefore the physical meaningful 



x/9« 2 -48rt+36c 2 +64 
6 

critical point (u cr i t ,v cr i t ,w cr it) is stable and is a late-time stationary attractor solution 



ranges are n < § ; c> and ~ 3 " 6 +] 4 " 2 +^ c 16 16 < Lo dm < |(3c + 1) - v»n--4»n+ 3 t> c -+b4 and in thig r£mge the 
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The dimensionlcss parameters u,v,w are drawn in figure 7 and we see that these are decreasing in progression 
of time, whereas u, w keep positive orientation and v keeps negative sign in late time evolution. Also v against 
u, w against v and w against u are drawn in figures 8, 9 and 10 respectively. Here v keeps negative sign when 
u and w increase. The 3D figure of it, v, w are drawn in figures 11 and 12 to get the visualization of the natures 
of it, v and w for 4 different initial conditions and a single initial condition respectively. The deceleration 
parameter q from figure 13 shows that q decreases from .5 to —1 which interprets the early deceleration and 
the late time acceleration generated by the VMCG in Kaluza-Klein Cosmology. 



V 



Fig. 11 




Fig.12 Fig. 13 



Fig. 11 and Fig. 12 show the 3D variations of the dimensionless parameters u,v and w for four different initial 
conditions and single initial condition. Fig. 13 shows the variation of deceleration parameter q against the time t. 



IV. CONCLUSIONS 

In this work, we have assumed the 5D Kaluza-Klein Cosmology for anisotropic universe where the universe 
is filled with variable modified Chaplygin gas (VMCG). The function B is assumed as a function of the scale 
factors a and b with some suitable choice B = B Q (a 3 6) ™. If the energy density and the pressure corresponding 
to the normal scalar field <\> and the self interacting potential V{<j>) can be associated with the energy density 
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and the pressure of the VMCG, the corresponding <f> and V have been calculated in terms of the scale factors 
a and b. The variations of V against a and b are drawn in figures 1 and 2 and have seen that V decreases 
as universe expands. Also from figure 3, we see that <p always increases as time goes on. Also we graphically 
analyzed the geometrical parameters named statefinder parameters in anisotropic Kaluza-Klein model. Fig. 4 
and Fig. 5 show the variations of r against a and b and the variations of s against a and b and Fig. 6 shows the 
variation of r against s respectively for some particular choices of the constants. From the figures, we conclude 
that r increases and s decreases as a and b increase. From figure 6, we see that s increases from some positive 
value to +oo as r decreases from finite positive range. Also after certain stage, s increases from -co to as r 
decreases from some positive to negative value. 

Next, we consider a Kaluza-Klein model of interacting VMCG with dark matter in the framework of 
Einstein gravity. Here we construct the three dimensional autonomous dynamical system of equations for this 
interacting model with the assumption that the dark energy and the dark matter are interact between them 
and for that we also choose the interaction term Q in the form AbHp, which is of course decreasing with the 
expansion of the universe. We convert that interaction terms to its dimensionless form and perform stability 
analysis and solve them numerically. The critical point (u cr it, v cr it, w cr it) have been found for the autonomous 
system. We obtain a stable scaling solution of the equations in Kaluza-Klein model and graphically represent 
solutions and have shown that the critical point is a late time stable attractor. The dimensionless parameters 
u, v, w are drawn in figure 7 and we see that these are decreasing in progression of time, whereas u, w keep 
positive orientation and v keeps negative sign in late time evolution. Also v against u, w against v and w 
against u are drawn in figures 8, 9 and 10 respectively. Here v keeps negative sign when u and w increase. The 
3D figure of u,v,w are drawn in figures 11 and 12 to get the visualization of the natures of u,v and w for 4 
different initial conditions and a single initial condition respectively. The deceleration parameter q from figure 
13 shows that q decreases from .5 to —1 which interprets the early deceleration and the late time acceleration 
generated by the VMCG in Kaluza-Klein Cosmology. 
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